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Abstract 

We adopt the so-called occupation number representation, originally used in quan- 
tum mechanics and recently considered in the description of stock markets, in the 
analysis of the dynamics of love relations. We start with a simple model, involving 
two actors (Alice and Bob): in the linear case we obtain periodic dynamics, whereas 
in the nonlinear regime either periodic or quasiperiodic solutions are found. Then 
we extend the model to a love triangle involving Alice, Bob and a third actress, 
Carla. Interesting features appear, and in particular we find analytical conditions 
for the linear model of love triangle to have periodic or quasiperiodic solutions. 
Numerical solutions are exhibited in the nonlinear case. 
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1 Introduction and preliminaries 



In a series of recent papers one of us (F.B.) used the framework of quantum mechanics, 
operator algebra and, in particular, of the so-called occupation number representation to 
discuss some toy models of stock markets [H EJ [31 H] . The main motivation for such an 
approach was that during the time evolution the main variables of a closed stock market 
can take only discrete values. This feature is well described by using the eigenvalues of 
some operators which describe these variables, the so-called observables of the market. 
Moreover, the closed stock markets we have considered admit conserved quantities, like 
the total number of shares or the total amount of cash, and these conserved quantities are 
well described in our framework. Here we want to show how the same general approach 
could be used in dealing with a completely different problem, i.e., the analysis of a love 
affair. In fact, it is natural to measure the mutual affection of the actors of our model using 
natural numbers (the higher the number, the stronger the love), and to think that some 
conserved quantities do exist in the game. It might be worth recalling that sophisticated 
mathematical tools have been used several times in the analysis of this problem, producing 
many interesting results which can be found in [51 El U\ IS 13 [10] , as well as in an extensive 
monograph, [IT] . For instance, in [9], some simple dynamical models involving coupled 
ordinary differential equations and describing the time variation of the love or hate in 
a romantic relationship are considered. In particular, a linear model for two individuals 
is discussed, and the extension to a love triangle, with nonlinearities causing chaotic 
dynamics, is also taken into account. 

It should also be considered that in the last few years a growing interest in classical 
application of quantum ideas has appeared in the literature. This involves very different 
fields like economics, [T2l |T3] , biology, [II] . and sociology, [151 [16] an d references therein, 
just to cite a few, and is a strong encouragement to carry on our analysis. 

The paper is organized as follows. In Section [21 we consider a first simple model 
involving two lovers, Alice and Bob, and we analyze the dynamics of their relationship 
starting from very natural technical assumptions. Both a linear and a nonlinear model are 
considered; then, the equations of motion are solved analytically (for the linear model) and 
numerically (for the nonlinear model), under suitable (and fairly good) approximations. 

In Section El we consider a model in which Bob has two relationships at the same 
time, and again we carry on our dynamical analysis. Also in this case we find an explicit 
solution for the linear model, while the nonlinear one is discussed numerically. Section [4] 
contains our conclusions, while, to keep the paper self-contained, the Appendix reviews 
few useful facts in quantum mechanics and occupation number representation. 
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2 A first model 



The first model we have in mind consists of a couple of lovers, Bob and Alice, which 
mutually interact exhibiting a certain interest for each other. Of course, there are several 
degrees of possible interest, and to a given Bob's interest for Alice (LoA, level of attraction) 
there corresponds a related reaction (i.e., a different LoA) of Alice for Bob. Now, let us see 
how this mechanism could be described in terms of creation and annihilation operators. 

Extending what has been done in [U [21 S] , we now introduce a± and a 2 , two inde- 
pendent bosonic operators. This means that they obey the commutation rules 

[ai,4] = aia{ - a\di = 18 i>k , (2.1) 

while all the other commutators are trivial: [a*, a^] = [aj, a k ] = 0, for all i and k. Further, 
let ip^ be the vacuum of aj, aj ip^ = 0, j = 1, 2. By using cp^ and the operators aj, we 
may construct the following vectors: 

:= -4=f ^ ¥W« := pff ® *>g, (2.2) 

where rij = 0, 1, 2, . . ., and j = 1,2. Let us also define Nj = aj aj, j = 1,2, and N = 
Ni + N 2 . Hence (see the Appendix), for all j, 

Nj aj ( Pm,n2 V^j -Q ^jVwi,n2) 

iVj a] v? nii „ 2 = (rij + 1) a)^ni,n 2 - 

Then, we also have N(p ni ^ n2 = (ni + n 2 )(p ni , n2 . As usual, the Hilbert space H in which 
the operators live is obtained by taking the closure of the linear span of all these vectors, 
for iij > 0, j — 1,2. A state over the system is a normalized linear functional w nijn2 
labeled by two quantum numbers n\ and n 2 such that oj nun2 (x) = ((p ni)Tl2 , x ip ni)n , 2 ) , where 
(., .) is the scalar product in % and x is an arbitrary operator on %. 

In this paper we associate the (integer) eigenvalue n\ of N\ to the LoA that Bob 
experiences for Alice: the higher the value of n\ the more Bob desires Alice. For instance, 
if n\ = 0, Bob just does not care about Alice. We use n 2 , the eigenvalue of N 2 , to label 
the attraction of Alice for Bob. A well known (surely simplified) law of attraction stated 
in our language says that if rii increases then n 2 decreases and viceversa0. Following the 
same general strategy as in [H 121 [31 H], this suggests to use the following self-adjoint 
operator to describe the dynamics (see Appendix) of the relationship: 

H = A (af 1 4 Ah + h.c.) . (2.3) 



lr This law has inspired many (not only) Italian love songs over the years! 
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Here, h.c. stands for hermitian conjugate, and M\ and M 2 give a measure of the kind of 
mutual reaction between Bob and Alice: if M\ is large compared to M 2 , Bob will change 
his status very fast compared with Alice. The opposite change is expected for M 2 much 
larger than M 1; while, for Mi close to M 2 , they will react essentially with the same speed. 
These claims will be justified in the rest of the paper. However, as it is already clear at 
this stage, it is enough to introduce a single index M rather than Mi and M 2 , since M 
plays the role of a relative behavior. For this reason we will choose the hamiltonian 



where A is the interaction parameter (which could also be seen as a time scaling parame- 
ter). The physical meaning of H can be deduced considering the action of, say, af 1 a\ on 
the vector describing the system at time t — 0, fni,n 2 - This means that, at t — 0, Bob is 
in the state ni, i.e. ni is Bob's LoA, while Alice is in the state n 2 . However, because of 
the definition of *p ni ,n 2 , a f f a \ V^i,^ which is different from zero only if M < m, is pro- 
portional to (f ni -M,n 2 +i- Hence, Bob's interest for Alice decreases of M units while Alice's 
interest for Bob increases of 1 unit. Of course, the Hamiltonian ( 12.4ft also contains the 
opposite behavior. Indeed, because of the presence of a 2 a{ in H, if n 2 > 1 we see that 
°»2 a { <fini,n 2 is proportional to v? ni+ A/,n 2 -i : hence, Bob's interest is increasing (of M units) 
while Alice looses interest in Bob. It is not hard to check that I(t) := Ni(t) + M N 2 (t) is 
a constant of motion: I(t) = 1(0) = Ni(0) + M iV 2 (0), for all tel. This is a consequence 
of the following commutation result: [H, I] = 0. Therefore, during the time evolution, a 
certain global attraction is preserved and can only be exchanged between Alice and Bob: 
notice that this reproduces our original point of view on the love relation between Alice 
and Bob. 

Now we have all the ingredients to derive the equations of motion for our model. 
These are found by assuming, as we have already done implicitly, the same Heisenberg- 
like dynamics which works perfectly for quantum systems, and which is natural in the 
present operatorial settings. Of course, this is a strong assumption and should be checked 
a posteriori, finding the dynamical behavior deduced in this way and showing that this 
gives reasonable results. More explicitly, we are assuming that the time evolution X(t) 
of a given observable X of the system is given by X(t) = e iHt Xe~ iHt . The equations of 
motion for the number operators Ni(t) and iV 2 (£), which are needed to deduce the rules 
of the attraction, turn out to be the following: 




(2.4) 



Ni(t) = i\M [al{t){ai{t)) 
N 2 (t) = -t\(ai(t)(ai(t)Y 



M 



M 



a 2 (t)(ai(tY) M ) 
a 2 (t)(ai(ty) M ) . 



(2.5) 
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By using this system, it is straightforward to check directly that I(t) does not depend 
on time. However, this system is not closed, so that it may be more convenient to replace 
(12.51) by the differential system for the annihilation operators ai(t) and a 2 (t): 



di(i) = -2AMa 2 (t)(a 1 (f)t) M " 1 
a 2 (t) — — i A (ai(t)) 



(2.6) 



Then, we may use the solutions of this system to construct Nj(t) = aj(t) cij(t), j = 1,2. 
Equations (12. 6p . together with their adjoint, produce a closed system. Of course, there 
exists a simple situation for which the system (12. 6 p can be exactly solved quite easily: 
M = 1. In this case, which corresponds to the assumption that Alice and Bob react with 
the same speed, system (12.61) is already closed and the solution is easily found: 

ai(t) = a\ cos(At) — ia 2 sin(A£), a 2 {t) = a 2 cos(At) — ia\ sin(Af). (2.7) 

Now, if we assume that at t — Bob and Alice are respectively in the ni'th and n 2 'th 
LoA's, the state of the system at t = is U} ni ,n 2 ( see Appendix). Therefore, calling 
rij(t) := L0 ni>n2 (Nj(t)), j = 1,2, we find that 

ni(t) = rii cos 2 (At) + n 2 sin 2 (At), n 2 {t) = n 2 cos 2 (At) + rii sin 2 (At), (2.8) 

so that, in particular, u nitTl2 (I(t)) = ni + n 2 , as expected. The conclusion is quite simple 
and close to our view of how the law of the attraction works: the infatuations of Alice and 
Bob oscillate in such a way that when Bob's LoA increases, that of Alice decreases and 
viceversa, with a period which is directly related to the value of the interaction parameter 
A. In particular, as is natural, setting A = in equation (I2.8p . implies that both Alice and 
Bob stay in their initial LoA's. The solution in (12. 8 p justifies, in a sense, our approach: 
the law of the attraction between Alice and Bob is all contained in a single operator, the 
hamiltonian of the model, whose explicit expression can be easily deduced using rather 
general arguments. The related dynamics is exactly the one we expected to find. Hence, 
the use of the Heisenberg equations of motion seems to be justified, at least for this 
simple model. A similar result was found for the linear model also in [9], under suitable 
assumptions about the parameters. 

Much harder is the situation when M > 1. In this case we do need to consider 
the adjoint of (12. 6p to close the system and, nevertheless, an exact solution can not be 
obtained. However, it is possible to generate a numerical scheme to find solutions of our 
problem, and this is the content of the next subsection. 
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2.1 Numerical results for M > 1 

The first remark is that, as stated above, the Hilbert space of our theory, H, is infinite- 
dimensional. This means that both Bob and Alice may experience infinite different LoA's. 
This makes the situation rather hard from a computational point of view and, furthermore, 
looks like a useless difficulty As a matter of fact, it is enough to assume that Bob 
(respectively, Alice) may pass through L\ (respectively, L 2 ) different LoA's (L\ and L 2 
fixed positive integers), which efficiently describe their mutual attraction. This means that 
our effective Hilbert space, "H e // ; is finite-dimensional and is generated by the orthonormal 
basis 

T = {(fin ■= Pn u n 2 , Tlj = 0, 1 Lj , j = 1 , 2} = U E K,} , (2.9) 

with obvious notation. Hence, dim(1-L e ff) = (Lj + 1)(L 2 + 1) is exactly the cardinality of 
/C. 

Calling H-eff the identity operator over % e ff the closure relation for J 7 looks like 
Yln&K Iv 9 " >< Pal = ^-eff- It 1S a standard exercise in quantum mechanics to check that 

(ai)k,k> ■=< Pk,a>i(Pk' >= \Zk[S ku k[-iSk 2 ,k' 2 , ^ ^ 

( a 2)k,k' ■ = < ¥k,a>2<Pk' >= V^2^2,fei-l ^fcl.fci) 

so that, by using the resolution of the identity in "H e //, we can produce a projected 
expression of the operators dj clS follows: 



(2.11) 



In other words, while the a/s act on "H, the related matrices A/s act on the finite dimen- 
sional Hilbert space "H e //- This means that, while the cij(£)'s are unbounded operators 
(which could be represented as infinite matrices), the A,(t)'s are (Li + 1)(L 2 + 1) x (Li + 
1)(L 2 + 1) matrices whose elements, for t = 0, can be deduced by (12. lip . Thus, these are 
really bounded operators and, as we will see, they all have the following property: for all 
choices of Lj a certain power Rj exists such that = 0, (A^) Rj = 0, condition which is 
not shared by the original cij. This fact is related to the approximation which produces 
7-teff out of the original "H. Now, system (12 .6p can be written in a formally identical way 
as 



a\ — > Ai 


== E 


< <Pk, 




\(fk >< (Pk'\ 


= E> 








k,k'eK 








k'eK 






a 2 ->■ A 2 


== E 


< <fk, 


a 2 <Pk> > 


\<fk >< tfk'l 


= E- 


/^2 






k,k'£K 








fe'e-ff 







i 1 (t) = -iAMA 2 (t)(A 1 (t)t)^-i 
i 2 (*) = -iA(Ai(t)) 
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where A\(t) and A 2 (t) are bounded matrices rather than unbounded operators. It should 
be stressed that the change a,- — > Aj does not destroy the existence of an integral of 
motion, which is clearly I e ff{t) = A/i(t) + MJ\f 2 (t), where A/}(t) = A^A^t). 

The numerical scheme is quite simple: we just have to fix the dimensionality of Heff, 
that is L\ and L 2 (i.e., the number of LoA's), the value of M in (12. 41) . and the matrices 
Aj(0). Then, by choosing a reliable scheme for integrating numerically a set of ordinary 
differential equations, we may construct a solution for the equations (12.1 2ft in a prescribed 
time interval. Fron now on, we take L\ = L 2 = K (K suitable positive integer), whereupon 
dim('H e ff) = (K + l) 2 , and we consider the orthonormal basis of T-i e ff i n the following 
order (the order is important to fix the form of the matrix!): 



Then we get 



{<A)0, Voi, • • • , Pok, Vio, <fiu, 



( o 



AM 



I 



4,(0) 



V 



K 


Ok 




Ik 


Ok 




Ok 




K ■ 


Ok 


K 


Vst K ■ 


Ok 


K 




X 


Ok 


• •• K 


K 


X 


••• Ok 


Ok 


K ■■ 


X 



K 


Ok \ 


Ok 


K 


Ok 


Ok 


Ok 


Ok 



VK - 11 K 



* / 



where Ok and \k are the null and the identity matrices of order K respectively, whereas 

/ ... \ 



x 





1 
y/2 0... 











\ o o ... y/ir^i o i 

It is easy to verify that ^(O)^- 1 ^ 0^ and ^(O)^" 1 ^ K , while Ai(0) K = A 2 (0) K = 
Ok- This simply means that, since because of our original approximation % — > H e ff, there 
are only K + 1 different levels; then, if we try to act more than K + 1 times on a certain 
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state, the only effect we get is just to annihilate that state. In other words, we can not 
move Bob or Alice to a — 1 or K + 1 LoA by acting with Aj or At since these states do 
not exist! However, this argument does not apply to the operators a,j in TL, since in the 
original Hilbert space there is no upper bound to the value of the LoA's of the two lovers. 

Suppose now that for t = the system is in the state u;„, where n = (ni,n 2 ). This 
means, as usual, that n describes the LoA of Alice and Bob. As we have done explicitly 
for M = 1, if we want to know how Bob's LoA varies with time, we have to compute 

in terms of the matrix elements. Analogously, to compute how Alice's LoA varies with 
time, we need to compute 

ra 2 (*) :=< ^n,A/" 2 (tK >= \\A 2 {t)^j 2 = i(^oy*)i 2 - 

The system (12. 12[) is numerically integrated for different choices of K by taking M = 2 
and A = 0.1. We use a variable order Adams-Bashforth-Moulton PECE solver [SU] as 
implemented in MATLAB®'s odell3 routine. 

In Figure [TJ the time evolutions of Alice's and Bob's LoA's in the case K = 3 with 
equal initial conditions n±(0) = n 2 (0) = 1 are displayed. Two clear oscillations in opposite 
phase can be observed: Alice and Bob react simultaneously but in different directions. 
This is in agreement with our naive point of view of Alice-Bob's love relationship, and 
confirms that the Heisenberg equations of motion can really be used to describe the 
dynamics of this classical system, also in this nonlinear (and likely more realistic) case. 
The horizontal dotted line on the top in Figure [I] (and also in Figures |2] and [3] below) 
represents the integral of motion I e ff(t); in some sense, it provides a check of the quality 
of the numerical solution. 

In Figure El the evolution of Alice's and Bob's LoA's with initial conditions (2,2) for 
different dimensions of "H e // is plotted. As it can be observed, when K = 3, the values 
of LoA's go beyond the maximum admissible value {K); this suggests to improve the 
approximation, i.e. to increase the dimension of 7i e ff- By taking K = 6, we find that 
Alice's and Bob's LoA's assume values within the bounds. Moreover, by performing the 
numerical integration with higher values of K [K = 7 and K = 8), the same dynamics 
as when K = 6 is recovered. We also see that, as in Figured! n\(t) increases when n 2 (t) 
decreases, and viceversa. 

The same main features are found in Figure El where the initial condition (1,3) and 
different values of K (K = 5, 6, 7, 8) are considered. The choice K = 5 turns out to be a 
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Bob 

Alice 

3.5 r 




time 

Figure 1: K = 3, M = 2: Alice's and Bob's LoA's vs. time with initial condition (1,1). A periodic 
behavior is observed. 

poor approximation; in fact, the values of LoA's can not be described in its T-L e ff- On the 
contrary, already for K = 6, the related % e ff turns out to be a fairly good substitute of "H, 
since the values of LoA's do not exceed the limits during the time evolution. Furthermore, 
the dynamics is stable for increasing values of K. 

Let us summarize the situation: as a consequence of our approximation, "H —> Heff, 
we need to replace I(t) with I e ff(t) = A/i(t) + MA/^t). As already noticed, because 
of (12.121) . I e ff(t) = 0, i.e., I e ff{t) is a constant of motion for the approximate model. 
Hence, it happens that if A/2 (t) decreases during its time evolution, A/i(t) must increase 
since I e ff(t) has to stay constant. In this way, for some value of t, it may happen that 
Ai(t) > K (see Figure [2] for K = 3, and Figure [3] for K — 5). This problem is cured 
simply by fixing higher values of K, as numerically shown in Figures |2] and El where the 
choice K = 6 allows us to capture the right dynamics, that remains unchanged also when 
K = 7 and K = 8. In a sense, the dimension of % e ff can be a priori fixed looking at 
the value of the integral of motion. In fact, if for some t, A/^t) = 0, then A/i(t) is equal 
to l e ff(0), which in turn depends on the initial conditions. Therefore, the dimension of 
% e ff must be greater than or equal to J e //(0); last, but not least, our numerical tests 
allow us to conjecture that, increasing the dimension of Heffi the resulting dynamics is 
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Bob 

Alice 




- Bob Bob 

Alice — - Alice 




Figure 2: Alice's and Bob's LoA's vs. time with initial conditions (2,2) and M = 2: K = 3 (top left), 
K = 6 (top right), K — 7 (bottom left), K = 8 (bottom right): K = 6 is already a good approximation. 
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Figure 3: Alice's and Bob's LoA's vs. time with initial condition (1,3) and M = 2: K = 5 (top 
left), K = 6 (top right), K = 7 (bottom left), K = 8 (bottom right): again, K — 6 is already a good 
approximation. 
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0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



frequency frequency 

Figure 4: Power spectra vs. frequency of numerical solutions for Alice's LoA (left) and Bob's LoA 
(right) for initial conditions (1, 1) (top: periodic dynamics) and (1,3) (bottom: quasiperiodic dynamics). 



not affected. 

It may be worth stressing that this feature only appears because of the numerical 
approach we are adopting here, and nothing has to do with our general framework. This is 
clear, for instance, considering the solution of the linear situation where no approximation 
is needed and no effect like this is observed. 

In this nonlinear case, the numerical results seem to show that a periodic motion 
is not recovered for all initial conditions, contrarily to what analytically proved in the 
case M = 1. In Figures |2] and El rather than a simple periodic time evolution (like that 
displayed in Figured]), a quasiperiodic behavior seems to emerge, as the plots of the power 
spectra of the time series representing the numerical solutions (see Figure HI bottom) also 
suggest. However, what emerges from our numerical tests is that, once the dimension of 
H e ff is sufficient to capture the dynamics, only periodic or quasiperiodic solutions are 
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obtained. This does not exclude that, for initial conditions requiring dimensions of the 
Hilbert space higher than those considered here, a richer dynamics (say, chaotic) may 
arise. Further numerical investigations in this direction are in progress. 



3 A generalization: a love triangle 

In this section we will generalize our previous model by inserting a third ingredient. Just 
to fix the ideas, we will assume that this third ingredient, Carla, is Bob's lover, and we 
will use the same technique to describe the interactions among the three. The situation 
can be summarized as follows: 

1. Bob can interact with both Alice and Carla, but Alice (respectively, Carla) does not 
suspect of Carla's (respectively, Alice's) role in Bob's life; 

2. if Bob's LoA for Alice increases then Alice's LoA for Bob decreases and viceversa; 

3. analogously, if Bob's LoA for Carla increases then Carla's LoA for Bob decreases 
and viceversa; 

4. if Bob's LoA for Alice increases then his LoA for Carla decreases (not necessarily 
by the same amount) and viceversa. 

In order to simplify the computations, we assume for the moment that the action 
and the reaction of the lowers have the same strength. Therefore, repeating our previous 
considerations, we take at first M — 1. The Hamiltonian of the system is a simple 
generalization of that in f!2.4[) with M — 1: 

H = A12 [a\ 2 a 2 + a i2 a\ J + A i3 (a\ 3 a 3 + a 13 a 3 j + A x (a{ 2 a i3 + a i2 c4 3 ) . (3.1) 

Here the indices 1, 2 and 3 stand for Bob, Alice and Carla, respectively, the A Q 's are real co- 
efficients measuring the relative interaction strengths and the different a a , a = 12, 13, 2, 3, 
are bosonic operators, such that [a a ,at] = 5 a ^1, while all the other commutators are 
zero. The three terms in H are respectively related to points 2., 3. and 4. of the 
above list. As usual, we also introduce some number operators: iV 12 = a| 2 a i2; describing 
Bob's LoA for Alice, iV 13 = a\ 3 ai 3 , describing Bob's LoA for Carla, N 2 = a\ a 2 , describ- 
ing Alice's LoA for Bob and N 3 = a\ a 3 , describing Carla's LoA for Bob. If we define 
J := N\2 + Nis + N2 + N3, which represents the total level of LoA of the triangle, this 
is a conserved quantity: J(t) = .7(0), since [H, J] = 0. It is also possible to check that 
[H, N u + N12] 7^ 0, so that the total Bob's LoA is not conserved during the time evolution. 
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The equations of motion for the variables a^'s can be deduced as usual and we find: 



i a 12 (t) = \ 12 a 2 (t) + \ t a 13 (t) , 
ia 13 (t) = X 13 a 3 (t) + Ai a 12 (t), 
ia 2 (t) = X 12 a 12 (t), 
k ia 3 (t) = \ 13 a n (t). 

This system can be explicitly solved and the solution can be written as 

A(t) = U' 1 exp (-iA d t) UA(0), 



(3.2) 



(3.3) 



where A is the (column) vector with components (ai 2 , a\ 3 , a 2 ,a 3 ), Ad is the diagonal matrix 
of the eigenvalues of the matrix 



A 



and U is the matrix which diagonalizes A. The eigenvalues of the matrix A are solutions 
of the characteristic polynomial 



/ 





Ai 


Al2 





\ 




Ai 








A13 






Al2 


















A 13 








I 



A 4 - (A? + \\ 2 + A? 3 )A 2 + A? 2 A? 3 = 0, 



(3.4) 



i.e., 



A ±)± = ±\l (A? + \\ 2 + A? 3 ) ± A /(A? + \\ 2 + X\ 3 f - 4A? 2 A? 3 . 
Since it is trivially 

-1, (3.6) 



A_ 



.± 



we have that the commensurability of all the eigenvalues is guaranteed if and only if 



with p and q nonvanishing positive integers, that is, if and only if the condition 

A? + A? 2 + A? 3 p , q 



(3.7) 



A12A 



12^13 



+ 

q p 



(3. 



holds true. 

Therefore, by computing Ni 2 (t), Ni 3 (t), N 2 (t) and N 3 (t), the solutions are in general 
quasiperiodic with two periods, and become periodic if the condition (13 .8p is satisfied. 
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Bob vs Carla 

Carla 




In Figure El we plot the solutions in the periodic case for L\ = L 2 = 3, ni 2 (0) = 0, 
^13(0) = 3, n 2 (0) = 713(0) = 2: Bob is strongly attracted by Carla, while both Carla and 
Alice experience the same LoA with respect to Bob. The parameters involved looks like 
these: A 12 = jq, A i3 = | and Ai = Jj, which satisfies (13. 8p for p = 2 and q — 1. 

The periodic behavior is clearly evident in both these plots. On the contrary, in 
Figure [6j we plot the solutions in the quasiperiodic case for the same L\, L 2 , ni 2 (0), 
n 13 (0), n 2 (0) and n 3 (0) as above. The values of the parameters are A i2 = ^, A13 = | and 
Ai ~ 0.0889, which satisfy condition (I3.8P for p = a/5 and q = 1. 

3.1 Another generalization 

A natural way to extend our previous Hamiltonian consists, in view of what we have 
done in Subsection 12. 1[ in introducing two different parameters M a , a = 12, 13, which are 
able to describe the different (relative) reactions in the two interactions Alice-Bob and 
Carla-Bob. We also assume that Bob is not very interested in choosing Carla rather than 
Alice, as long as one of the two is attracted by him. For these reasons, the Hamiltonian 
looks now like 

H = A 12 (V 12 ) Ml2 a 2 + a^ 12 4) + A13 ((« f i 3 ) Ml3 % + «i 3 12 4) + *i («i2 «is + «i2 « f 13 ) ■ 
Also in this case an integral of motion does exist, and looks like 

J := N 12 + N 13 + M 12 A^ 2 + Mi 3 iV 3 , 
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Bob vs Alice 

Alice 



Bob vs Carla 

Carla 




Figure 6: (Xi = L 2 = 3: Alice's and Bob's LoA's vs. time with initial condition (0, 2) (left) and Carla's 
and Bob's LoA's vs. time with initial condition (2,3) (right). Quasiperiodic behaviours are observed. 



which reduces to J if M 12 
in (HOD, i.e., 



M 13 = 1. The equations of motion also extend those given 



% d l2 (t) = A 12 Mia (a\ 2 (t)) M ^ l a 2 (t) + X 1 a 13 (t), 

ia lz {t) = A 13 Mia (alaW)** 13 " 1 a sW + A i M*), 
za 2 (f) = Ai 2 a 12 (t), 
ia 3 (t) = Ai 3 ai 3 (t). 



(3.9) 



System (13.91) is nonlinear and can not be solved analytically unless if M 12 = M 13 = 1, 
as shown before. Therefore, we adopt the same approach we have considered in our first 
model. 

In particular, we take L\ = L 2 = 2, consider the orthonormal basis J 7 = {<fijki}, 
where the indices i,j,k,l run over the values 0,1,2 in such a way the sequence of four- 
digit numbers ijkl (in base-3 numeral system) is sorted in ascending order, whereupon 
dim('H e ff) = (Li + 1) 2 (L 2 + 1) 2 = 81. In such a situation, the unknowns ai 2 (t), ai 3 (t), a 2 {t) 
and a 3 (f) in the system (13. 9p are square matrices of dimension 81, and their expression 
at t — are given by: 
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Figure 7: (L a = 2, M 12 = 1, M13 = 2 (a = 12,13,2,3): LoA vs. time of: Bob vs. Alice (continuous 
line), Bob vs. Carla (dashed line), Alice (dashed-dottcd line), Carla (dotted line) with initial condition 
(1, 0, 2, 2). Periodic behaviours are observed. 



where 
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Moreover, in order to consider a situation not very far from a linear one, we choose 
Mia = 1, Mi 3 = 2, Ax = 0.3, A12 = 0.2, A i3 = 0.007. 

The numerical integration has been performed by using once again the odell3 routine 
of MATLAB®, and the existence of the integral of motion has been used in order to check 
the accuracy of the solutions. 

The choice Ai 3 = 0.007 has an immediate consequence, which is clearly displayed by 
the plots in figures ITlfTOl The time evolution of N 3 (t) is very slow and essentially trivial 
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2.4 
2.2 




Figure 8: (L a = 2, M 12 = 1, M13 = 2 (a = 12,13,2,3): LoA vs. time of: Bob vs. Alice (continuous 
line), Bob vs. Carla (dashed line), Alice (dashed-dottcd line), Carla (dotted line) with initial condition 
(1, 2, 1, 2). Periodic behaviours are observed. 



in the time interval where the system is numerically integrated: Carla keeps her initial 
status. This does not imply, however, that Bob does not change his LoA concerning 
Carla, as the figures ITIfTOl also clearly show. Another interesting feature, which is again 
due to the nonlinear aspect of the dynamics (and it was not present in the linear model), 
is related to the fact that, for certain time intervals, Bob's LoA for both Alice and Carla 
can increase. This is, in a certain sense, unexpected: the Hamiltonian does not present 
any explicit ingredient responsible for this feature, which, in our opinion, appears just 
because of the nonlinearity of the interaction. Finally, we also notice that the figures 17101)1 
display a periodic behavior, even if no analytic proof of the existence of this periodicity 
was produced, up to now. 

It is worth noticing that, contrarily to what happened in Section 2 (for certain initial 
conditions the dimension of H e ff needs to be increased in order to capture the right 
dynamics), all the numerical results obtained here suggest that, once the dimension of 
~H e ff is fixed (and this is done when the initial conditions are given), the time evolution 
of the various LoA's can be described inside that 'He//; &.e., there is no need to enlarge 
the dimension of the effective Hilbert space. This could be related to our choice of the 
parameters, which, as already stated, produces a system with small nonlinear interactions. 
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Figure 9: (L a = 2, M 12 = 1, M13 = 2 (a = 12,13,2,3): LoA vs. time of: Bob vs. Alice (continuous 
line), Bob vs. Carla (dashed line), Alice (dashed-dottcd line), Carla (dotted line) with initial condition 
(2, 0, 0, 2). Periodic behaviours are observed. 

For this reason, we do not expect this to be a general feature of the model: on the contrary, 
we expect that, changing significantly the values of the parameters, the need for increasing 
the dimensionality of Heff wm arise also here. Further analysis in this direction is in 
progress. 

4 Conclusions 

In this paper we have shown how to use quantum mechanical tools in the analysis of 
dynamical systems which model a two- and a three-actor love relationship. Depending 
on the parameters which describe the system, linear or nonlinear differential equations 
are recovered, and exact or numerical solutions are found. These solutions show that a 
nontrivial dynamical behavior can be obtained, and a necessary and sufficient condition 
for the periodicity or quasiperiodicity of the solution in the linear love triangle model has 
been found, (13. 8p . 

It is worth stressing that, despite the quantum framework adopted in this paper, all 
the observables relevant for the description of the models considered here commute among 
them and, for this reason, can be measured simultaneously and no Heisenberg uncertainty 
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time 



Figure 10: (L a = 2, M 12 = 1, M 13 = 2 (a = 12, 13,2,3): LoA vs. time of: Bob vs. Alice (continuous 
line), Bob vs. Carla (dashed line), Alice (dashed-dottcd line), Carla (dotted line) with initial condition 
(2, 1, 0, 2). Periodic behaviours are observed. 

principle should be invoked for these operators. 

A final comment is related to the oscillatory behavior recovered in this paper for the 
different LoA's. Indeed, this is expected, due to the hamiltonian approach considered 
here. In a paper which is now in preparation, [17] . one of us is considering the possibility 
of describing, using the same techniques, other processes in which some time decay may 
occur because of possible interactions of the lovers with the environment: this, we believe, 
could be a reasonable method to construct a more realistic model, where also different 
psychological mechanisms can be taken into account. 

Appendix: A few results on the occupation number 
representation 

We discuss here a few important facts in quantum mechanics and second quantization, 
paying not much attention to mathematical problems arising from the fact that the op- 
erators involved are quite often unbounded. More details can be found, for instance, in 

[mug. 
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Let % be an Hilbert space, and B(H) the set of all the bounded operators on %. 
Let S be our physical system, and A the set of all the operators useful for a complete 
description of S, which includes the observables of S. For simplicity, it is convenient 
to assume that A coincides with B(Ji) itself, even if this is not always possible. This 
aspect, related to the importance of some unbounded operators within our scheme, will 
not be considered here, even because the Hilbert space H e ff considered in this paper is 
finite-dimensional, which implies that the operators in A are bounded matrices. The 
description of the time evolution of S is related to a self-adjoint operator H = W which 
is called the Hamiltonian of S, and which in standard quantum mechanics represents the 
energy of S. We will adopt here the so-called Heisenberg representation, in which the 
time evolution of an observable X e A is given by 

X(t) = e iHt Xe- iHt , (A.l) 

or, equivalently, by the solution of the differential equation 

dX ^ ie im [H,X]e- iHt = i[H,X(t)}, (A.2) 



dt 

where [A, B] := AB — BA is the commutator between A and B. The time evolution 
defined in this way is usually a one-parameter group of automorphisms of A. 

An operator Z G A is a constant of motion if it commutes with H . Indeed, in this 
case, equation (1A.2|) implies that Z(t) = 0, so that Z(t) = Z for all t. 

In our paper a special role is played by the so-called canonical commutation relations 
(CCR): we say that a set of operators {a£, a\,£ = 1,2, satisfy the CCR if the 

conditions 

[a e , a*] =6 en l, [a t , a n ] = [a{, o) n ] = (A. 3) 

hold true for all £, n = 1,2, ...,L. Here, 1 is the identity operator. These operators, 
which are widely analyzed in any textbook in quantum mechanics (see, for instance, |18j ) 
are those which are used to describe L different modes of bosons. From these operators 
we can construct hi = a\ai and iV = Y2e=i which are both self-adjoint. In particular, 
hi is the number operator for the £-th mode, while is the number operator of S. 

The Hilbert space of our system is constructed as follows: we introduce the vacuum 
of the theory, that is a vector </?o which is annihilated by all the operators ae: aiipo = 
for all i = 1, 2, . . . , L. Then we act on ip with the operators a\ and their powers: 

VW..,,* := . , \ M r {a\r ■ ■ ■ (4)" L ^o, (A.4) 
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rii = 0, 1,2, . . . for all /. These vectors form an orthonormal set and are eigenstates of 
both hi and N: n£(p ni ,n 2 ,...,n L = nt<p nim ,...,n L and N^ ni) n 2 ,...,n L = N(p nun2t ^ nij , where 
N = J2e=i n i- Moreover, using the CCR we deduce that ni{Q>i<Pni,n2 l ... ) n L ) = ( n e ~ 

m,n 2 ,...,n L ) an d fl£ \ 0'£ ( Pni,n2 l ...,nL 

= (n e + l)(aj 

sons, the following interpretation is given: if the L different modes of bosons of S are 
described by the vector ip ni ,n 2 , -,n L , this implies that ri\ bosons are in the first mode, ri2 in 
the second mode, and so on. The operator he acts on (p ni ,n 2 ,...,n L and returns ne, which is 
exactly the number of bosons in the I— th mode. The operator N counts the total number 
of bosons. Moreover, the operator ae destroys a boson in the ^-th mode, while a\ creates 
a boson in the same mode. This is why ae and a\ are usually called the annihilation and 
the creation operators. 

The Hilbert space H is obtained by taking the closure of the linear span of all these 
vectors. 



The vector (f ni ,n 2 ,...,n L in (1A.4j) defines a vector (or number) state over the algebra A 

as 

w 7li ,712,. ni,n2,—,TiLi -^(Pni,ri2 hl 

>, (A.5) 

where ( , ) is the scalar product in the Hilbert space "H. As we have discussed in [H |2l [31 H] , 
these states are used to project from quantum to classical dynamics and to fix the initial 
conditions of the system considered. 
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